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I — I Summary 

^ This study deals with spectral analysis of fluid-loaded vibrating structure. It was recently observed [8] 
I . in a numerical study on a high order perturbation method under heavy fluid loading that a loaded 
c/2 ! vibrating plate results, not only in the classical frequency shift of the in vacuo single resonance (in 

a 

Q ' a result of the loading, a single in vacuo resonance frequency of the structure is transformed into 



both the real part because of the fluid added mass and the imaginary part because of energy lost 
by radiation into the fluid), but also in an increase in the number of the resonance frequencies : as 



> 



o 



a multiple resonance frequency. Here we show that this phenomenon is said to be a refinement of 
the Sanchez's classical result [T^ (paragraph 9.3) where it was established, using asymptotic analysis, 
P i that in the case of a light loading conditions "the scattering frequencies of a fluid loaded elastic 
structure (ie the resonance frequencies) are nearly the real eigenfrequencies of the elastic body alone 
and the complex scattering frequencies of the fluid with a rigid solid" . 

A theoretical explanation of the observed phenomenon of multiple resonance is given using classi- 
cs! cal results of the distribution of zeros of entire functions. It is established that every single in vacuo 
resonance frequency of a rectangular plate is transformed into an infinite number of resonances when 
^ the fluid-loading is accounted for. 
O PACS no. 43.20.Ks, 43.20.Tb, 43.40.Rj 

o 

> 1 Introduction 

X 

^ _ It is well known that the need to accurately describe fluid structure interactions increases the com- 
■ - - ' putational cost of the various numerical methods used. Most numerical methods can be used to deal 
with the resolution of linear systems of simultaneous equations. When the loading conditions are 
taken into account, the size of the matrices involved increases considerably, and they become full 
and frequency dependent. Any method of reducing these drawbacks is therefore most welcome; and 
one of the best methods available is that based on the perturbation approach. Asymptotic analysis 
can be conduced when the loading is light, as in the case of a metal structure surrounded by air. 
This approach consists in introducing a small parameter e, which is the ratio between the surface 
mass density and the fluid density. Using perturbation expansions [S], one can then construct an 
approximate solution, based on the in vacuo eigenmodes, which is not only easier to calculate than 
the solution of the exact problem but also leads to a better understanding of the various phenomena 
involved [3]. 

The light fluid loading perturbation method (involving high density structures in contact with 
a small density fluid) was recently extended [8] to cases where the perturbation parameter becomes 
large (which corresponds to a light and/or thin structure in contact with a high density fluid). 
During this study multiple resonance processes were observed where a single resonance mode can 
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have several resonance frequencies. The spectrum of the operator can no longer be described in terms 
of the discrete spectrum based on a countable series of resonance mode/resonance frequency pairs. 
This typically nonlinear phenomenon is closely linked to the concept of non-linear frequency modes 
(0)' where the non linearity depends on the time parameters (such as those involved in porous 
materials) rather than on the more usual geometrical parameters. This frequency non-linearity is 
related to problems which are written in the form 0{uj)U = S, where u is the unknown factor and S 
is the source term. 0{u}) is an operator which depends non linearly on the angular frequency uj, as 
in the case of porous environments or in the context of vibroacoustics, where the coupling depends 
non-linearly on the frequency via the Helmholtz equation Green's kernel. 

Here it is established that this behavior is a general property of at least fluid-loaded simply 
supported rectangular plates. The main results obtained in our previous study [8jJ are recalled in 
paragraph [21 in particular, the method for calculating the spectrum of a fluid loaded vibrating 
structure based on a high order perturbation expansion is recalled and the resonance frequency 
equation to be solved is given. A theoretical description of multiple resonance in the case of a 
simply supported plate loaded with a non vanishing density fluid is given in paragraph [21 It is then 
established, using classical results on the connection between the order of an entire function and the 
distribution of its zeros, that any resonance mode in this structure will have an infinite number of 
resonance frequencies. Numerical examples of a steel plate in contact with air or water are given 
paragraph [H The conclusions and the possible extensions of this work are presented in paragraph [51 



2 Statement of the problem 

Consider a finite elastic structure occupying a domain S, the behavior of which is described by 
a differential operator A (in the case treated in this communication it is a simple bi-Laplacian 
operator). This structure is loaded with a perfect fluid at rest extending to infinity. The equation of 
the fluid-loaded structure is given by the classical integrodifferential equation 

AU{M, uj) - pshuj^ (u{M, ^) - e ^ U{M', uj)G{M; M', uj)dM'^ = F{M) (1) 

where e = pf/psh is a small parameter in the case of a metal structure (volume mass density ps and 
thickness h) in contact with air (volume mass density p/). G{M] M', uj) is the Green's function of the 
Neumann problem in the Helmholtz equation outside the surface S. F{M) is the source. To facilitate 
reading, the weak formulation of the problem is introduced. One defines the radiation impedance 
P{U,V,uj) = Jj. Jj.U{M')G{M; M' ,uj)V*{M)dMdM', where U{M) and V{M) are functions defined 
on the surface S, a{U, V*) stands for the potential elastic energy of the structure, {U, V*) is the usual 
inner product and pshLv'^{U, V*) is the kinetic energy. Then the weak formulation reads: find U{M) 
such that for every V{M) one has 

a{U, V*) - AM {{U, V*) - ePiU, V*, u)) = {F, V*). (2) 



2.1 Resonance modes and eigenmodes for a fluid- loaded structure 

The eigenmodes Um{M,u) and eigenvalues Am{uj) (or equivalently the eigenpulsations Umi^) given 
by Am{uj) = pshoj^{uj)) are the non zero solutions of the homogeneous weak formulation that is 
find Um{M,uj) and uj^iuj) such that for every V{M) 

a{U^, V*) - pshuiioo) [{Um, V*) - e(3{Um, V*, u)) = 0, (3) 

/?([/, V, cu), that characterizes the energy loss by radiation, depends on the frequency. Then the 
eigenmodes, well suited to deal with harmonic regime, are complex and depend on the frequency. 
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One introduces the resonance modes, that allow to compute the transient regime in a natural 
way; they are complex and do not depend on the frequency [3]. The resonance modes Um{M) and 
resonances pulsations Um, given by = Pshu"^, are the non zero solutions of the following problem: 
find Um{M) and Um such that for every V{M) 

a{lJ^, V*) - pshul ((f/„, V*) - el3{iJm, V\ u^)) = 0. (4) 

These two kind of modes, that are identical when the loading is not took into account, have a 
strong link as each resonance pulsation is the solution of the resonance equation 

'^m(^^m) = C^^, (5) 

and each resonance mode is an eigenmode calculated at the corresponding resonance pulsation. 

U^{M) = U^{M,u^). (6) 

It is worth noting that in the eigenmodes are much more easier to compute than the resonance 
modes. 



2.2 High order perturbation method 

When the loading parameter e is small, the eigenmodes and eigenvalues can be calculated using a 
perturbation method. The following perturbation expansion is introduced into the weak formula- 
tion ([3D 

Um{M) = u(^)(M) + eU!:^\M) + --- + e'Ut\M) + --- 

where U^\M) and A^^) are the eigenmodes and eigenvalues of the elastic structure in vacuo. To 
facilitate the reading, one poses: (3"^"'{uj) = [3 (u^\ U^'>*,u?j. It can then be easily estabhshed that 
one has 



to the order 1 in e:A(^)(cj) = A(^)/3""^(cu), 



to the order 2 in e:A(^)(^) = A^ + ^ " .(o) /^"^"(^) 

\ p^m Ap 

and so on to the higher orders. Next, as previously shown [8], by grouping the terms and by a 
re-summation technique, one shows that the high order perturbation expansion is given by 

A(o) 

which remains valid even for higher values of e. It is worth noting that, compared to the classical 
perturbation expansion of the eigenvalues given by A„i{uj) ^ A(^) (l + e/5'^'"(cj)), the high-order 
approximation ([7]) is particularly interesting as it gives very precise results for the same numerical 
cost. 

The calculation of the eigenmodes and eigenvalues does not requires any particular effort and, 
for a given u, each eigenpulsation is the solution of the following eigenpulsation equation 

Equation shows that the calculation of the various eigenpulsations do not requires any com- 
putational cost but the evaluation of multiple integrals /J™"* (a;). 
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2.3 The resonance equation 

By introducing in equation ([5]) the high order expansion of the eigenpulsation given by equation ([8]), 
the resonance equation ([5]) is approximated as 



This shows that the calculation of the resonance pulsations is completely different than that of 
the eigenpulsations since (a;) has a rather complicated dependence on the frequency and solving 
equation iQ requires the search for the zeros of non-convex complex function which is a challenging 
task. A numerical example given in the paper [8J has shown a curious phenomenon of multiple 
resonance. It has been made a comparison of the evolution of the resonance frequency in the mode 
(1 — 3) of clamped steel plate loaded by water when thickness vary using an exact solution and 
the approximated one (by the high order perturbation expansion). That particular example has 
showed that the resonance frequency is not unique and at least two resonance frequencies could exist 
simultaneously for a given geometrical configuration. 

To have a better understanding of this non-usual phenomenon the roots of the resonance equa- 
tion ([9]) are now studied in a more general manner using classical results in the framework of the 
distribution of zero of entire functions [HI E]. In the following the resonance equation ([9]) is written as 



where z = ujm and zom = ^m- The zeros of this function are the resonance pulsations that are looked 
for. 

3 Multiple resonances 
3.1 Methodology 

Now all the problem focuses on the existence and the counting of the roots of J"*™ (2;) = (1 — e/5™'™(2;)) — 
Zq^. Let us presents the methodology used in the rest of paper. More or less, the theory of the zeros 
of entire functions (analytic functions in the whole complex plane) is based on the Great Picard's 
Theorem [11] which states that if /(z) has an essential singularity at a point then on any open 
set containing fi^z) takes on all possible values, with at most one exception, infinitely often. The 
exception is needed as showed by the function exp(l/2;) that has an essential singularity at z = 
but never attains as a value. The second point is that every entire function is either a polynomial 
or has an essential singularity at infinity. Now if one can prove that the resonance function f"^"^{z) 
is an entire function, that is not a polynomial and that the value is not the exception needed by 
the Picard's theorem, then the resonance function f"^"^{z) has an infinite number of zeros. Then for 
every mode, an infinite number of resonance exists. 

Now to characterize complex function f{z) as an entire function, the simplest way is to study 
its development as a power series of the form f{z) = YI^q' Cnz'^. When lim„_^oo a/c^ = 0, f{z) is 
an entire function. To show that it is not a polynomial, one can study its order. Let us recall that 
if f{z) is an entire function, the function M(r) = max|^|=r \ f{z)\ will increase indefinitely together 
with r. The rate of growth of M(r) can be characterized by comparing it with function exp(r). 
By definition, the order (of growth) p is the smallest number such that M(r) < expr''"'"'' for every 
> and for r > roi^r]). From this, it can be easily seen that the order is given by the formula 
p = limsup^^j^ In In M(r)/ In r ; as examples, a polynomial in z is of order 0, and the exponential 
function is of order 1. The order p of an entire function depends on the asymptotic behavior of the 
coefficients of its series expansion. One has 1/p = lim„^oo iiif(ln(l/|c„|)/(nlnn)) (^HJ); or in other 
words, if Cn oc l/(n!)° then p = 1/a. 



c:;^(l-6/?— (^J)^^W2. 



(9) 





(10) 
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One can simplify again the problem under consideration. The resonance function can be written 
as f'^"^{z) = — Zqj^ — ez'^ j3"^"^ (z) , as a sum of a polynomial {z^ — Zq^) of order and a function 
which is the product of a polynomial {ez"^) also of order and the radiation impedance (3^™'{z) which 
order is to be determined. The distribution of zeros of f^^{z) can be deduced from that of 13"^"^ {z) as 
long as e has a non zero value. If (3"^"^{z) is a polynomial, the resonance function is also a polynomial 
with a finite number of roots. If I3^^{z) is an entire function (in the case considered below it is of 
order 1) with essential singularity at infinity, then f^^{z) is also an entire function with essential 
singularity at infinity and by the Picard's theorem f"^"^(^z) can takes on all possible values at infinity, 
with at most one exception, infinitely often. 

The main results of this paper is, as shown below on the particular example of a simply supported 
plate, that the radiation impedance and the resonance function are entire functions of exponential 
type (that is of order one and normal type [5]) that both have an infinite number of roots. 

3.2 The resonance function of a simply supported baffled plate as an 
entire function of the frequency 

Now let us consider a thin simply supported rectangular plate (dimensions a x b with thickness 
h) made of isotropic material with Young's modulus E, Poisson's coefficient u and density ps- In a 
vacuum, the eigenmodes are given by U^l {x,y) = 2/ \/ab sin rmrx / a sin mry / b and the eigenpulsations 
by ul^l = 7T^JD/{psh) [{m/af + {n/bf] with D = Ehyi2{l - z/^). When the plate is baffled, the 



radiation impedance P{U^l,Uj^}^* ,uj) is given as [lOj 

~ / ~ \ db P^ P^ €^~^ 

f^i^mni^mn 1^) — I I I i siu mTTX siu WTT?/ siu mTTx' siu nvT?/' — —dxdydx'dy', {11] 

TT JO JO JO JO d 



where d = \J a?[x — x'Y + h^iy — y'Y and c is the speed of waves in the fluid. In the following, the 
pulsation dependence in uj is changed into the equivalent wavenumber dependence in k = uj / c. Since 
the in vacuo modes involve a double index, the notation for the radiation impedance is simplifled as 
/3(t/i°), U!^l\ u) = /?'""(A;). The resonance function is denoted by /"^"(A;) =^{1- e/?""(A;)) - k^^^, 
with komn = ^^^l/c. 

It is worth noting that /5™'"(/c) and f"^'^{k) decrease in the upper complex plane (with '^{k) > 0, 
the term exp{tkd) has an exponential decrease when \k\ — > oo) and conversely grow in the lower 
complex plane, then all the roots of these two functions are located in the lower complex half plane. 
As the time dependency had been chosen as exp(— zcut), this ensures that the movement of the plate 
remains flnite over time [5]. 

By making a change of variable [7], it can be easily shown that I3^^{k) is given by 

^-(A;) = / Gm{X)G^{Y)--dXdY, (12) 

TT JO JO U 



where D = \J d^X"^ + b'^Y'^ . The functions Gm{p^^ and Gn{y) are given by single analytic integrals 
Gmi^) = Jq sin rmT(x + x') sin mnx'dx' with G^ix) = (sin rrnr x + mn(l — x) cos rmrx) /{2mn). By 
making a change of variables from rectangular to polar [1], with < ^^q, = arctan(6/a) < tt/2, the 
radiation impedance f|T2l) is written as 

pmn^,^^ = _1 j^^^^^k, e)de - - r'^ J2„„(A;, e)de, (13) 



7r JQ TT 
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with 



Jimn{k,9) = ('°" H^n{R.e)e'^''dR, (14) 

b 

J2mn{k,e) = f H^niR,e)e''''dR, (15) 

JO 

HUR,e) = gJ^)gJ^]. (16) 



By expanding the exponential as a power series, the integrals ( IT4l) and (fT5|) are given as 

Jlmn{k,0) = 22 m \ Hmn{R-, d)R dR, 

V. Jo 



e=oo 



J2mn{k,0) = ^ '^~m~ I Hmn{R,0)R^dR 



as Hmn{R, 0)R^ is a continuous integrable function of R g]0, a/ cos6'[ and R g]0, 6/ sin6'[, the mean 
value theorem ensures that there exist Ri g]0, a/ cos6'[ and R2 g]0, h/ sin6'[ such that 



/ HmniR, 9)R^dR = HmniRl, 0)R{, 
Jo 

/■ills 

/ Hmn{R, S)R dR = HjnniR2, G)R2- 

Jo 



10 

Similarly as HmniR,^^) is a continuous integrable function of 9 G]0,7r/2[, there exist 9i g]0,6'q,[ and 
92 ^]9a,TT/'^[ such that 



I Jlmn{k,9)d9 = ^ ——Hmn{Rl,9l)R{, 

Jo ^^Q I. 

I J2mn{k,9)d9 = ^ Hmn{R2,92)R2- 



One then obtains: 



r"(A:) = Ur.n{R.. 9,) E + ^^^.(^2, 92) ^ , (17) 

f3"^^{k) is then an entire function of the variable k. 
3.2.1 Order of the radiation impedance 

To calculate the order of the impedance radiation, one studies the comportment of the function 
Mpm7i[r) = max|fc|=r |/5™"(A;)| for r —>■ oo. Let us begins with successive integrations by parts of the 
integrals given by equations f|T^ and f|T5|) . Introducing the results given in the appendix leads to 

JiUk, ^) = + (;^^2™n(^^) - ^ {hl^^{9)e^'^ - hl^^{9)) 

1 d'^Hmn{R,()) tkR 



+W^o dR^ ' 

1 d^Hmn{R,&) ikR 
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As one studies the comportment of max|fc|^oo \P"^^(yk)\, in the equations (fTS!) and (fT9l) . the first two 
power terms in k are neghgible compared to the exponential ones when \k\ — >■ oo; similarly, only the 
first exponential term has a contribution when \k\ oo. Then 

1 



max 



-^ik- 



>d9 



\k\=r 

Now by expanding the exponentials in equation fl20|) as power series, one obtains 



(20) 



max — — 

\k\=r \k\^ 



l=oo 
f=0 




COS 9 



sin 9 



d9 



[21) 



As h\^^^{9) and a/ cos 6* are continuous integrable functions for 9 G [0, ^q] and h\^^^^{9) and 6/ sin^ 
are continuous integrable functions for 9 G [^Q,7r/2]. One applies the mean value theorem for the 
two integrals that appear in equation fl2T|) : there exist two values 9i and 92 such that 



/ a 



Vcos 

f— 

\ sin^ 



d9 



d9 



^3mn(^2) 



COS^l 

sin 9o 



,0<9i<9a 



,9a<92< n/2. 



(22) 
(23) 



Let us denote ci = a/ cos 6*1 and C2 = b/ sin 6*2- It is easy to see a < ci < \/a?~+lP and 
6 < C2 < \/a^ + W. Mpmn[r) is then written as 

M^mn (r) 



max 


1 


^3m?i( 


^=oo 
£=0 


(i/cci)^ 


^3mn( 


£=oo 
92)E 


{lkC2Y 


\k\^ 


£! 




max 

|fc|=r 


1 


^3mn( 


ei)e*^"^ 


+ ^3mn( 


^2)e*^"^ 







(24) 

The previous result shows that the function Mi3mn[r), characterizing the growth of the radiation 
impedance, is the sum of two exponential functions divided by a polynomial of the fourth degree. 
Its order is easy to obtain from classical results on the sum and product of entire functions (see for 
example [11], chapter VII paragraph 6). There is two possibilities for characterizing the growth of 
Mjjmn{r) for r sufficiently large: the first is that if a denotes a finite positive number and C a positive 
constant one has Mpm.n{r) Cexp(crr) which implies that the radiation impedance is of order 1, 
the second is that Mpmn^r) = which implies that the radiation impedance is of order 0. From the 
expressions of the functions hl^^{9) and hl^j^{9) given in the appendix, it is obvious that for a 7^ 6 
and for m 7^ n one cannot have ^smnl^i)^*'^'^^ 



'^Imni^^)^''^'^^ ■ If, for example Ci > C2, then 



max 



he 



ifcci 



|fc|=r 

If a = 6 and m = n care must be taken; in this case ^q, = 7r/4 and a little algebra shows that 



Then the maximum of f3"^'^{k) is 



cos 9 



d9 



Ma 



(r) 



max 



k\=r \k\- 



sin 9 



d9. 



(25) 



(26) 



(27) 



Then /3"^^{k) is an entire function of order 1. The resonance function /"^"^{k) = k'^{l — ef3"^^{k)) 
k1 is an also an entire function of order 1 for all e 7^ 0. 
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3.2.2 Distribution of zeros of the radiation impedance and of the resonance function 

In addition to the order, the relations fl24l) . (!25|l and (!27|) give a more refined information on the 
distribution of zeros of the radiation impedance. An entire function that grows as exp{ar), with 
< cr < oo, is an entire function of order one and normal type a. Such functions are known as 
entire function of exponential type (EFET) [B] . The resonance function and the radiation impedance 
are EFET. They both possess an additional property on the distribution of their zeros. They are 
functions that belong to the Cartwright class C [6j. A EFET f{z) belongs to the class C if 



— at < oo, 

1 + ^2 



(2^ 



\w\ 



max(0, log or equivalent ly as log'' 



\w\ 



where the function log \w\ is defined as log 
l/(27r) /o^^ log \w - exp{te)\de. 

A function f(t) bounded on the real axis satisfies log"*" < \f(t)\. To prove that f{z) belongs 
to the class C it is sufficient to prove that is bounded for every t E IR and belongs to class C. 

The radiation impedance given by equation (fT3|) can be majored, for k E IR, by: 



|/3™"(A;)| 



< 



4 
vr 
4 



tt/2 



b 

sin 6 



H^Jr,e)e''''dr\de 



\H^n{r, e)\dr]de+ i'^ (j^'" \H^^{r, 9)\dr\ do) . (29) 



As Hmn{r,d) is continuous integrable function of r and 9 the two integrals that appear in the right 
member of equation (l29l) are finite. Then (3'^^{k) is bounded for every k real and belongs to class C. 

Showing that the resonance function /"""(fc) = k'^{l — e/5''""(fc)) + fegmn belongs to class C requires 
a more refined estimation as, because of the k'^ term, it is unbounded at infinity. Nevertheless, as 
/3"^^{k) is bounded for every k, the resonance function is also bounded for every finite real value k, 
then for every finite value A 

-A 1 + t2 

It remains to estimate /f^ i°g+^l/J(*)l rft. As seen in the previous paragraph, the radiation impedance 
for real values of ^ oo has the following asymptotic approximation 



-dt < oo. 



(30) 



/3™"(fc) 



[tky 



(31) 



The two integrals that appear in this equation : 



h{k) 



^LnW^^de 



can both be approximated for k ±oo by the stationary phase method [T]. One obtains for k — ±oo 

3mn (0)\ 



hik) 



hi 



2ka 



2a3 V 2ka 
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Now the radiation impedance equation has the following approximation for k ±00 




For k ±00, /?™"(A;) tends to zero and the resonance function can be approximated by /™"(A;) ^ k"^. 
For k — > ±00 one has log^ |/™"(^)| = log |A;^|. Then for every value of A 

/ ^ / ^' rft ^ 2 / -^rft = 2 / ^^rft = 2 ^ 7 + Oil/A") < 00. (33) 

Combining results of equations (|30ll and (l33l) shows that < 00. Then both 

impedance radiation and resonance relation belong to the Cartwright class C. 

For < a < TT let us denote by n+(r, a) the number of zeros of the function f"^"'{z) in the sector 
{z : 1^1 < r, I arg(z)| < a} and n_{r,a) the number of zeros in {z : \z\ < r, |7r — arg(z)| < a}. By 
the theorem of Cartwright and Levinson (theorem 1, page 127 of |6]) the density of the set of zeros 
n±{r, a)/r has the limit \imr^oon±{r, a)/r = a/n. This result is a very strong one, since not only it 
ensures that the resonance function possesses an infinite number of zeros but also that "almost all" 
zeros lie in the neighborhood of the real axis. 

It is worth noting that the zeros are counted with their multiplicity order. Then the zeros of 
the resonance function appears as conjugate pairs and then must be counted twice. Actually, the 
resonance function /'^"'(fc) = A;^(l — e/5'""(A;)) + ^omn can be approached for \k\ 00 by /"-"(A;) = 
fc^(l — e2/i3^.^(6'i)//cV"'i + ^omn; ^^^u if is a solution of /"""-(fc) = 0, it is easy to see that —k*, the 
opposite of its conjugate satisfies f"^"'{—k*) = 0. 



4 Numerical examples 

In this section, some numerical examples of the multiple resonance phenomenon are presented. One 
considers a steel [E = 200 Gpa, u = 0.3 and ps = 7800 Kg/m^) simply supported rectangular plate 
in contact with air (with density pf = 1.2 Kg/m^ and sound wave celerity c = 340 m/s) or water 
(with density pf = 1000 Kg/m^ and sound wave celerity c = 1500 m/s). The dimensions of the plate 
are indicated in the chosen examples. 

As it is impossible to obtain a closed formula for Jimn{k, 6) and J2mn{k, 6) given by equations 
and ( |T5l) . a numerically efficient approximation of these functions has been constructed. As shown 
below, this approximation preserves the order properties of the involved functions. Since one always 
has < 6a < vr/2, the functions to be integrated are expanded with respect to the angular variable 
9 to the order 1. This gives 

,, ^, za^fc^ + 2m^7r^ — zafcm^TT^ — (— l)™2e*"'^m^7r^ , , 

Jimn{k,e) ^ jim{k) = a — , 34 

4[ak — mTTY[ak + rmiY 

ih^k^ + 271^7^"^ -ihkn'^Ti'^ - i-lY2e^^n^n'^ , , 

J2mnik,e) ^J2nik) = b — . 35 

4[bk — m[Y[nk + m^y 



After reporting the two expansions fl34|) and fl34j) in the radiation impedance given in equation f[T3|) . 
the following approximation of (3"^^{k), denoted by /3™'"(fc), is obtained 

/3""(A;) = --UiUk)0a+j2n{k){n/2 - 6^)). (36) 
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It is worth noting that the accuracy of the approximation of i3"^'^{k) decreases as the mode order 
increases. However, at low orders (say m < 3, n < 3) it remains sufficient. It can easily be see that 

■ fl ^ / \ ±3t + 17171 

jim{k = ±mTC a) = a — — , (37) 

lomTT 



jlmik = 0) 



1 - (-1)™ 



±3z + nvr 

j2n[k = ±1171 h) = b—- , (3^ 

loriTT 



J2n{k = 0) = b 



i-(-ir ^ 



Then as jim{k) and j2n{k) are entire functions of k, $"^^{k) is also an entire function of k. It still 
remains to estimate the order of /3™"(A;). For given m and n and for \k\ — >■ cxo, one has 



_(_l)m^2^ie 
^IfrajUk) = 

lim j2n{k) 



and 



.m 



2 




hm r^'ik) = 27r ( (-1)™ ^„e"^^ + (-1)"^ ( ^ - ) e^"' ) TT- (39) 



It is worth noting that the approximation of the radiation impedance given in equation (!36|) or in 
equation (1391) presents a similar comportment than the exact one: a sum of exponential functions 
divided by a polynomial of the fourth degree. To confirm that the approximation keeps the order 
of (3"^'^{k) to one, classical results on the sum and product of entire functions are used again. First 
it is worth noting that the entire function (— l)™m^/a^^Qe*°'^ + (— l)"'n^/6^ (7r/2 — 9a) e*^^ cannot be 
identically zero for any value of a, b, m or n when \k\ —>■ oo. To see this more clearly, one can deal 
separately with the cases a = b and a^^ b. 

The case a = b. One has 6^ = 7r/4 and /3'""(fc) ^ ^ {{-ly^m^ + n^) e'"''^ becomes the 
product of the exponential function e*"'^ of order 1 and the inverse of a degree 4 polynomial in k, 
which is of order 0. $"^^{k) is therefore an entire function of order 1 and of type a. 



The case a ^ b. Let us take for example a > b. Then with \k\ — oo one has max 



> 



max 



; it is worth noting that this is true in the half plane '^{k) < 0, whereas in the half plane 
Q{k) > 0, /3"^"{k) tend rapidly to zero. In the half plane Q'(fc) < 0, one therefore has $"^^{k) ~ 
^ (^(— l)™'m^^Q,e*°'^^ ^ which is also the product of the exponential function e*"'^ of order 1 and the 

inverse of a degree 4 polynomial in k, which is of order 0. I3"^"'{k) is therefore an entire function of 
order 1 and of type a. It is worth noting that if 6 > a its type is b. In other words, one can assume 
that a = max(a, b). 

Figures ([H) and ([2]) give contour plots for the amplitude of the radiation impedance (3"^"'{k) and 
its app:oximation j3"^^{k) in the mode m = 1 = 1 in figure ([1]) and in the mode m = 1 n = 3 in 
figure ([2]). In these curves, the dark spots correspond to the zeros of the radiation impedances. It is 
worth noting that the lower the mode order, the better the approximation /3™"(/c) becomes. 

Both radiation impedance (3"^^{k) and resonance relation f"^"'{k) belong to the Cartwright class 
C and then all, but a finite number of, their zeros must lie in the neighborhood of the real axis. 
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Figure 2: Contour plots of the amplitude of the radiation impedance in the mode m = 1, n = 3 with 
a = 1 m and 6 = 0.7 m . Left: exact P^^{k), right: approximate P^^{k). 
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As shown above the type a of /3™"(A;) is given by a = max(a, b). Then the type of is given 

the type of A;^(l — eP"^"'{k)) + fcomn! for e 7^ it is given by cr = max(a, 6). By the Cartwright 
and Levinson theorem, the density of resonances limr^ooin+ir, a) /r + n^{r,a)/r) is then given by 
2max(a, b)/n. 

To verify this, the simplest way is to count the zeros of using the argument theorem on 

the disk C with radius K centered on the origin. It is clear from equations (157|) and (155]) that /3"^"'{k) 
and f"^"'{k) have no pole inside C. The number of its zeros inside this disk, in the case of the mode 
(m — n), is denoted by njmn{K). It is given by the integral njmn{K) = l/(2z7r) f''^^{k)/f'^'"'{k)dk 
with k = Ke^^ . This integral is computed numerically without difficulty. As shown before, the zeros 
occur in combined pairs of negative imaginary part, then each zero must be counted twice; then the 
density of the zeros given by {n{K,a)/K = n^{K,a)/K + n_{K,a)/K has the asymptotic value 
\imK^oon{K,a)/K = 2nj!mn{K)/K. 

For example, let us take air for the fluid and let the plate dimensions be a = 1 m, 6 = 0.7 m 
and thickness h = 0.002 m. It can easily be established that in the mode m = 1, n = 1, one has 
nj_^_^{0.3) = 2, ?2j^^(12) = 2 and nj^^(30) = 8. In the mode m = 1, n = 3, one has n^_^^{2) = 2, 
nj^g(lO) = 2, nj^^(30) = 8. As the radius K of the disk C becomes large one obtains for the density 
nfmn{K) /K the following values 

a = 1 m, 6 = 
a = 1 m, 6 = 
a = 1 m, 6 = 
a = 0.5 m, b 
a = 0.5 m, b 
a = 0.1 m, 6 
a = 1 m, 6 = 

This shows that lim^- 
are a validation of the theoretical results obtained in the paragraph [31 

The resonance relation of the fluid loaded plate has an infinite number of zeros. Therefore, while 
in a vacuum, each mode of a simply supported plate will have only one resonance frequency (and its 
opposite), under loading condition each mode will have an infinite number of resonance frequencies. 
This result is confirmed in figures [3] and H] where the dark spots that appear in the curves correspond 
to the zeros of the approximate resonance value relation fmn{k) = 0. The results presented in 
figures [5] and H] also show that except for the zero close to the real axis, the zeros of the resonance 
value relation corresponding to the plate loaded with water are very similar to those given by the 
radiation impedance. 

Now let us return to the initial question of the existence or non-existence of multiple resonances 
in fiuid-loaded plates. In every boundary condition, one has to look for the resonance pulsations Comn 
such that uJmni^ ~ ^/^<Jw(^)) " ^mn = 0- Although it has not been proved that such a comportment 
is true under various boundary conditions, it seems reasonable to assume that under boundary 
conditions of all kinds every resonance mode in a plate has an infinite number of resonance frequencies 



-- 0.7 m,K = 20; * ^; = 1/3.3333 
K 

-- 0.7 m, = 200; ^ ^; ^ = 1/3.125 
K 

-- 0.7 m,K = 20000; ^ ^; ^ = 1/3.14159 

K 

= 0.5 m,K = 20000; ^ ' = 1/6.28319 

K 

= 0.1 m, = 20000; * ^ = 1/6.28319 

K 

= 0.1m,K = 20000; ^ = 1/31.4159 

K 

rifuiK) nfisiK) nf32(K) 
-- 0.7 m,K = 20000; ^ ^ ' = ^ ^ ' = ^ ^ ' = 1 /3. 14159 

K K K 

^oon{K,a)/K = 2njmn{K)/K = 2 max(a, 6)/7r. All these numeric values 
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20 30 40 50 



Figure 3: Contour plots of the amplitude of the resonance value equation /c^(l — e/3'""(/c)) — k^^^ in 
the mode m = 1, n = 1 in the case with a = 1 m , 6 = 0.7 m and h = 2 mm. Left: steel plate in 
contact with air (e ~ 0.08), right: steel plate in contact with water (e ~ 64). 




Figure 4: Contour plots of the amplitude of the resonance value equation k^{l — 
the mode m = 1, n = 3 in the case with a = 1 m , b = 0.7 m and h = 2 mm. 
contact with air (e ~ 0.08), right: steel plate in contact with water (e ~ 64). 



e$^-{k)) - in 
Left: steel plate in 
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Figure 5: Contour plots of the amplitude of the resonance value equation k^{l — e/3™"(A;)) — k^^^ in 
the case e = 10^^ with a = 1 m , 6 = 0.7 m and h = 2 mm. Left: in the mode m = 1, n = 1, right: 
in the mode m = 1, n = 3. 

as e 7^ (the radiation impedance given by relation f[T^ seems to be an EFET for a large variety of 
functions Gm and 

As shown by Sanchez P2|, under very light loading condition (as in the case of a thick steel 
plate in contact with air, for example), there exist only the resonances of the structure alone (one 
in each mode shifted toward the complex plane) and the complex scattering resonances of the fluid 
in contact with the rigid solid. This contradictory situation is only an apparent one since for each 
resonance one has a root located very close to imaginary axis (that corresponds to the resonance 
of the in vacuo mode shifted through the complex plane by the coupling with the fluid with added 
mass and added damping) and an infinite series of roots corresponding to the zeros of the radiation 
impedance. For e — > 0, that is the case considered by Sanchez, all resonances sets introduced by 
the radiation impedance of each mode merge together and form a unique set of resonances for the 
resonance function. The demonstration seems very difficult but a numerical example of this merging 
is showed in the figure Oin the modes m = 1, n = 1 and m = 1 n = 3 in the arbitrary case e = 10~^. 
In this case, the structure has only two sets of resonances in accordance to the results by Sanchez, the 
first one is the set of resonances of the structure that are very close to the real axis plane (composed 
by the set of resonances of each modes of the structure with small imaginary part) and the second 
which is the resonance set of the rigid structure alone (corresponding to the merged zeros of the 
radiation impedance with large negative imaginary part). 

Now under very heavy loading condition (as in the case of a very thin plate in contact with heavy 
density fiuid ), what is the relevant relation ? This seems to be an open question. Although the 
resonance relation k'^(l — eP™'^{k)) — k^}^ = gives good approximation to the spectrum under heavy 
loading [S], it has to be supplemented by terms accounting for the cross modal impedance, which 
make it almost impossible to conduct the above analysis. 

5 Conclusions 

The aim of this paper was to show that the loading of a vibrating plate with a non- vanishing density 
fiuid transforms each resonance of the structure into an infinite number of resonances. This result 
opens up new prospects. In particular, it would be interesting to obtain similar estimates on the 
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vibratory modes occurring in plates with various boundary conditions, shape or loading (e.g. for non 
baffled plate) or in mechanical structures of other kinds, although these may not be easy to obtain. 
Another question which arises is how to determine under what loading conditions (moderate or 
heavy) a coupled system, such as a thin steel plate in contact with water, will have these properties. 
Another important aspect worth investigating is the validity of modal expansions in the case of 
resonance modes having two or more resonances frequencies such as those studied in this paper. It 
would also be interesting to develop experimental device which would make it possible to test the 
occurrence of the behavior described here. 
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A Derivatives of the function Hmn{r, 6) with respect to r 

The function Hmn{f, &) is given by 



1 \ ( rcos9\ frmrrcos9\ f rmrrcosO^ 

/ rsin6'\ /n7rrsin^\ fmrr sin theta\ 
nn [^--^) cos [—^) [ 1 j 

A tedious calculation (done without difficulty under any symbolic program) shows that the deriva- 
tives with respect to r of Hjnn{r, 0) up to order 3 taken at the points r = 0, r = a/ cos 9 and r = b/ sin 9 
are given by : 

1 

4 




r=0 



[Hmn{r,9)l^^ 

COS 6 

[Hmnir,9)l^^ 

sin 

dHmn{r,9y 
dr 

dr 

' dHmnir,9) 
dr 



= 

= 

= 

= 



d^H„,n{r,9y 
' d^Hmn{r,9y 

Qj.2 

d'Hmn{r,9y 



r=0 



. b 
'sine 



r=0 



d^H^n{r,9) 
dr^ 



TT^ (b'^m'^ cos^ 9 + a^n^ sin^ 9 









1 



' rn^ cos^ 9 ri^ sin^ 9 ' 



-71 



+ 



l)"'m^TT cos^ 9 
2a^bn 



63 



, _ / / an7rtan^\ . / anTrtan^' 
cos 9 IriTT cos I I + sm I 



/ amTta.n9\ 
-amr cos : sm 9 



{-iyn^7rsin^9 
2ab^m 



^ ( ( bmn cot \ / bmn cot 9 ' 
a sm 9 rmr cos + sm 

V V « y V « . 



—brmr cos ^ cos 



/ bmn cot 9 ' 



